SIGNAL AVERAGING USING GATING SIGNAL OBTAINED FROM 
AUTOCORRELATION OF INPUT SIGNALS 

Field of the Invention 

The present invention is directed to signal averaging to alleviate noise in signals having 
periodic or almost periodic waveforms, e.g., those taken in pulse volume measurement, and more 
particularly to such signal averaging in situations in which a gating signal is unavailable. 
Description of Related Art 

The measurement of peripheral blood flow is important in medicine, since there are many 
specific diseases in which peripheral blood flow is impaired, e.g., diabetes and atherosclerosis. 
Also, the peripheral blood flow is altered as the total cardiac output is increased or decreased. 
Cardiac output is particularly important in patients who are under anesthesia, are in the post- 
operative state, or are critically ill or unstable. As blood flow from the heart falls, the peripheral 
blood flow is dramatically reduced to preserve flow to the brain and vital organs. 

Blood flow to an extremity can be measured painstakingly and invasively by dissecting 
out the main blood vessels to the limb (e.g., brachial artery in an arm) and encircling it with an 
electromagnetic flow probe. That is clearly not a technique suitable for clinical use. It is 
therefore desired to measure peripheral blood flow non-invasively. 

Peripheral impedance (or conductance) plethysmography is a technique for non- 
invasively measuring peripheral blood flow by measuring peripheral pulse volume, which is the 
small change in the volume of a limb segment occurring within the cardiac cycle. The technique 
works by obtaining a raw pulse volume analog signal and applying a selective signal-averaging 
algorithm to the raw pulse volume signal. The technique is described in U.S. Patent No. 
4,548,211 to Marks. 



117622.00103/35583333v2 



1 



In the technique as currently practiced, the raw pulse volume analog signal is obtained by 
measuring the electrical impedance (or conductance) of a limb segment with an electrode such as 
that of Fig. 1. The electrode 102 is made of a flexible material 104, so that it can be wrapped 
around the limb segment. The flexible material is configured to define a connecting portion or 

5 member 106, which is insulated from direct electrical contact with the patient, and two members 
108 for being wrapped around or otherwise applied to the limb segment. Each of the two 
members 108 contains an outer current electrode 110 paired with an inner voltage electrode 112. 
An electrical connector 114 allows the outer current electrodes 110 and the inner voltage 
electrodes 112 to be connected to a source of current and a voltage-measuring device, 

10 respectively. 

An alternating current on the order of 1 ma amplitude and 40 kHz frequency is applied to 
the two outer current electrodes, while the inner voltage electrodes are used to measure the 
voltage resulting from the applied current. The ratio of the amplitude of the voltage waveform to 
the amplitude of the current waveform is the limb impedance, Z. Measurements of Z over time 
15 provide the baseline impedance Z 0 of the limb segment and the pulsatile change AZ of the 
impedance. Knowing the resistivity p of the blood and the distance L between the two inner 
voltage electrodes, the change in volume A V can be calculated as: 
AV=pL 2 AZ/Z 0 . 

The resistivity p is either calculated or approximated from the patient's hematocrit. An 
20 improvement to the above technique is taught in U.S. Patent Application No. 10/392,308, filed 
March 20, 2003, by Smith et al. 

The signal-to-noise ratio in computer-assisted plethysmography can be improved through 
selective signal averaging. The ECG signal is used as a gating signal. More specifically, the 
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gating signal is based on detection of a QRS complex, which is the series of deflections in an 
electrocardiogram that represent electrical activity generated by ventricular depolarization prior 
to contraction of the ventricles. When a QRS complex is detected, a finite time window is 
opened during which the analog input signal must remain within certain amplitude boundaries. 
5 If the analog input signal remains within the defined amplitude boundaries, the waveform is 
saved for averaging. If not, it is discarded. 

However, the above technique has the following drawbacks. To obtain a time reference 
for the averaging process, an ECG signal must be provided. That requires either additional leads 
to be applied to the patient or an ECG signal to be provided from another patient monitor. Either 
10 of those possibilities increases the number of leads to be applied to the patient and the 
complexity of the overall equipment to be used. 

In an entirely separate field of endeavor, the mathematical technique of autocorrelation 
has been developed to detect features of data such as periodicity and randomness. 
Autocorrelation can be considered as the correlation or convolution of a signal with itself, so that 
15 the signal is compared to a time-shifted version of itself. Mathematically, the autocorrelation of 
a signal of finite duration can be expressed as 

00 

R(r)= \f{t)f(t + r)dt. 

-co 

If fif) is real, the autocorrelation function R(r) is real and even. Thus, the above expression is 
equivalent to 

CO 

20 R(r)= \f{t)f{t-t)dt. 

-00 

The autocorrelation o f random noise takes the form o f a spike, since random noise i s 
similar only to the non-time-shifted version of itself. On the other hand, the autocorrelation of a 
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periodic signal, in which the time-shifted version goes into and out of phase relative to the non- 
time-shifted version as the time is shifted, is also periodic. Since a short signal is similar to itself 
only while it lasts, its autocorrelation drops to zero once the time shift exceeds the lifetime of the 
signal. 

5 The use of autocorrelation in digital signal processing is well known, as is the design of 

digital signal processors to carry it out. However, insofar as known, the technique of 
autocorrelation has not previously been applied to the above-noted limitations of 
plethysmography. 
Summary of the Invention 

10 It will be readily apparent to those skilled in the art to which the present invention 

pertains that a need exists in the art to eliminate the need for an ECG signal or other separately 
supplied gating signal. 

It is therefore an object of the present invention to provide a system and method for 

deriving the gating signal from the input signal itself. 
1 5 Another object of the invention is. to use such a derived gating signal to locate individual 

waveforms in the input signal. 

Still another object of the invention is to use such a derived gating signal to determine 

time periods in which the waveforms are tested to determine whether they remain in a required 

amplitude range for averaging. 
20 To achieve the foregoing and other objects, the present invention is directed to a system 

and method for obtaining a timing reference for the averaging process. It is based upon the 

autocorrelation (or equivalently, the convolution) signal processing technique. As an illustrative 

example, consider an analog signal that contains two or more successive plethysmographic 
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waveforms. If the noise content is not excessive, the peak of the autocorrelation function applied 
to that analog segment will occur at a time interval equal to the distance between the two peaks 
of the successive plethysmographic waveforms. 

In practice, a segment of analog input including more than two waveforms will be 
5 acquired and analyzed. The autocorrelation function does not contain a single, precise peak 
because there is natural variability of the distance between successive plethysmographic 
waveforms. This is due to cardiac R-R interval variability, which in turn is due, at least in part, to 
respiration. 

In that case, the autocorrelation function provides "first pass" information that determines 
10 the range over which the distance between the plethysmographic waveforms varies. Then, a 
single waveform can be isolated, reversed, and "slid over" the entire waveform segment using 
the convolution method (that determines how the single waveform should be positioned in time 
relative to other single plethysmographic waveforms) to facilitate the averaging process. This 
may be repeated for each single waveform. 
1 5 An additional b enefit i s that the amplitude o f the c orrelation function c an b e u sed for 

determining signal-to-noise quality and can provide an additional criterion for waveform 
acceptance for, or exclusion from, averaging. This technique may be applied to other signal 
averaging tasks in which a gating signal is either not available or difficult to obtain. 
Brief Description of the Drawin2S 
20 A preferred embodiment of the present invention is disclosed in detail with reference to 

the drawings, in which: 

Fig. 1 illustrates a plethysmographic electrode according to the prior art; 
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Fig. 2 is a block diagram of a system on which the preferred embodiment of the method 
of the invention is implemented; 

Fig. 3 depicts an analog plethysmographic signal; 

Fig. 4 shows the signal of Fig. 3 compared with a time-shifted version of itself for 

5 autocorrelation; 

Fig. 5 shows an isolated waveform of the plethysmographic signal; 
Fig. 6 shows the isolated waveform of Fig. 5 after time reversal; 

Fig. 7 shows the original signal compared with the reversed isolated waveform for 
convolution; and 

I o Fig. 8 is a flow chart of the operation of the preferred embodiment. 

Detailed Description of the Preferred E mbodiments 

A preferred embodiment of the present invention is set forth in detail hereinafter with 

reference to the drawings. 

Fig. 2 shows a block diagram of a device in which the preferred embodiment can be 
1 5 implemented. The device 200 differs from conventional plethysmographic devices and from that 
of the aforementioned Marks patent and Smith et al application in that, inter alia, it includes a 
digital signal processor for performing autocorrelation and convolution functions. 

The electrode 102 and the processing device 220 are connected by way of electrical 
connectors 114, 222. In the processing device 220, a current source 224 applies alternating 
20 current to the two outer current electrodes to induce a voltage in the two inner voltage electrodes. 
A voltage-measuring device 226 measures the induced voltage and supplies the measured value 
to a microprocessor 228 or other suitable processing device. The microprocessor 228 is in 
communication with a memory 230, which can be a ROM, an EEPROM, or other suitable non- 
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volatile memory, and which stores software routines to perform the operations to be described 
below. The microprocessor 228 is also in communication with a digital signal processor (DSP) 
chip 234 that performs DSP functions, such as convolution and autocorrelation. The 
microprocessor 228 and the DSP chip 234 perform the operations explained below to average the 
5 waveforms; then the microprocessor 228 calculates V and outputs the calculated value to any 
suitable output 232. 

The device 200 operates in a manner that will now be explained with reference to the 
waveforms of Figs. 3-7 and the flow chart of Fig. 8. An analog signal selected by the device 
(Fig. 8, step 802) has an appearance substantially as shown in Fig. 3. As shown, the analog 
1 0 signal 5(0 has a plurality of plethysmography waveforms defined by peaks P\,Pi,Pi, P*, Ps, Pe, 
P 7 , and troughs T u T 2 , h, T 4 , T 5 , T 6 . If 5(0 is autocorrelated with a time-shifted version of itself, 
S(t-x), as shown in Fig. 4 and specified in Fig. 8, step 804. Then, if the noise in 5(0 is not 
excessive, the resulting function of t has peaks occurring at a time interval equal to the time 
interval between successive peaks of S(t). Each waveform to be used in the averaging should be 
1 5 tested t o d etermine w hether i ts a mplitude r emains i n a r ange b etween a m inimum a cceptable 
amplitude A mi „ and a maximum acceptable amplitude A max . 

However, in practice, the time interval between successive peaks of 5(0 has some 
variation due to cardiac R-R interval variability, which in turn is due, at least in part, to 
respiration. Therefore, the autocorrelation function does not contain a single, precise peak. 
20 To overcome that difficulty, the autocorrelation function is not taken as the final answer, 

but instead provides first-pass information about the range over which the distance between the 
successive peaks of 5(0 varies. By use of that first-pass information, a single waveform W(t), t = 
0 through to, can be isolated, as shown in Fig. 5, as specified in Fig. 8, step 806, and as tested in 
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Fig. 8, steps 806 and 808, to determine whether its amplitude remains in the boundaries required 
for averaging; if not, it is discarded in Fig. 8, step 810. That waveform, if not discarded, can 
then be reversed to form W(t 0 -t), as shown in Fig. 6. The reversed waveform can be convolved 
with 5(0, or "slid over" S(t), as shown in Fig. 7 and specified in Fig. 8, step 812. The 
5 convolution of two waveforms shown in Fig. 7 can be mathematically expressed as 

C(r)= \S(t)W(r-t)dt. 

The peaks (local maxima) of the convolution as a function of x give the best matches between 
the peak of W(t) and the peaks of the other single plethysmographic waveforms included in 5(0- 
Thus, the sliding or convolution can be used to determine the time window during which the 
10 input signal must remain within the amplitude boundaries mentioned above, as determined in 
Fig. 8, step 814. The sliding or convolution can also be used to determine how the waveform 
Wit) should be positioned in time relative to the other single plethysmographic waveforms for 
the averaging operation, as specified in Fig. 8, step 816. The isolated waveform can be matched 
and averaged with multiple other waveforms in the signal, as specified in Fig. 8, step 818. Each 
15 single waveform in S{i) can be isolated, so that the sliding can be repeated for each single 
waveform, as specified in Fig. 8, step 820. The matching and averaging process ends in step 
822, whereupon an averaged signal is available for calculating V. 

The present invention is not limited to plethysmography, but instead can be applied to 
any signal averaging task in which it is difficult or impossible to obtain a gating signal. In 
20 particular, the present invention can be used with regard to any signal that is periodic or almost 
periodic and that has considerable noise. One example is analysis of a waveform indicating 
shaking in a device such an automobile. Since the same problems with signal averaging may 
arise with such a waveform as those noted above with a plethysmographic waveform, the present 
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invention is applicable in such situations as well. Of course, those skilled in the art who have 
reviewed the present disclosure will readily appreciate the other situations in which the present 
invention can be used. 

While a preferred embodiment of the present invention has been set forth above, those 
skilled in the art who have reviewed the present disclosure will readily appreciate that other 
embodiments can be realized within the scope of the invention. For example, as explained 
above, the present invention can be applied to other signal averaging tasks in which it is difficult 
or impossible to obtain a gating signal. Also, disclosures of specific waveforms are illustrative 
rather than limiting. Therefore, the present invention should be construed as limited only by the 
appended claims. 
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